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Abstract—Brightness temperature change of observation target in one imaging period will 
introduce image blur to the retrieved image of geostationary interferometric microwave sounder 
(GIMS), which uses rotating circular array to realize time-sharing sampling other than snap-shot to 
reduce system complexity but meanwhile resulting in a relatively long imaging period. In this paper, 
an interpolation method utilizing time sequence of visibilities collected by time-sharing sampling is 
presented. The method is verified through a simulation system that simulates GIMS’s observation 
process in viewing a series of near real case brightness temperature maps modeled by 
FNL/WRF/RTTOV method. Results show that the interpolation method can de-blur the image and 
reduce imaging error to the degree of snap-shot imaging. 


INTRODUCTION 


Nowadays, some geostationary earth orbit (GEO) microwave instruments at 50-56GHz and 183GHz have been 
proposed for its ability to retrieve atmospheric temperature and humidity profiles from radiometer observations 
and help investigate internal structure of tropical cyclone. Compared with low earth orbit sounding, sensors 
working at GEO orbit can continuously monitor the full earth disk, thus performing real-time observation for fast 
changing weather such as tropical cyclones. However, to achieve a moderate spatial resolution around oxygen 
band (53GHz) in GEO observation, a quite large real aperture antenna is needed. Furthermore, deploying and 
scanning the large traditional real aperture antenna poses a big challenge for GEO satellite application. To 
overcome the problem, the technique of interferometric aperture synthesis has been put forward. It uses signals 
intercepted by multiple small antennas to yield angular response characteristics of a much larger antenna, 
realizing high spatial resolution while relieving system complexity. With this idea, some instruments have been 
proposed including GeoSTAR by Jet Propulsion Laboratory, NASA [1]; GAS by European Space Research and 
Technology Center, ESA [2]; and geostationary interferometric microwave radiometer (GIMS) by National 
Space Science Center, Chinese Academy of Sciences [3]. Fig 1 shows an artistic view of GIMS instrument. 


BLURRING PROBLEM AND DEBLURRING METHOD 


As is shown in fig 1, GIMS uses a rotating circular array [4] and is supposed to be working in the time-sharing 
mode. Fig 2 displays the instantaneous sampling coverage in frequency domain at two different time and can 
help illustrate time-sharing operating mode. The rotating circular array system collects part of visibilities in one 
integral period thus obtaining the incomplete instantaneous sampling coverage in frequency domain at time t1 


shown in fig 2(a). As the circular array system rotates with time, assuming that it rotates an angle © from time 
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tl to t2, it will collect another set of visibilities at time t2, leading to the instantaneous sampling coverage at time 
t2 positioned at an angle © to the sampling coverage at time tl, which is shown in fig 2(b). Since the final 
image can’t be retrieved until complete sampling coverage is obtained in frequency domain, and samples in 
frequency domain are conjugate symmetric, the rotating circular array needs to rotate half a circle to complete 
collecting all the visibilities needed for image retrieval, therefore leading to a longer imaging period than that of 
the snap-shot operating mode, which collects all the visibilities needed for image retrieval in one integral period. 


Fig 1. Artistic view of GIMS instrument 
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Fig 2. Instantaneous sampling coverage at time t1 and time t2 ((a) time t1 (b) time t2) 


However, atmospheric parameters are under continuous variation, resulting in time-variant observational 
brightness temperatures of the full earth disk. As a consequence, brightness temperature variations of 
observational target in the relatively long imaging period might introduce some blurring to the retrieved image 
for GIMS. In addition, according to target brightness temperature modeling for GIMS in [5], brightness 
temperature variation is particularly obvious in tropical cyclone area, hence a method is expected to be employed 
to eliminate blurring and restore the retrieved image. 

An intuitive method is the interpolation method utilizing time sequence of visibilities collected by 
time-sharing sampling mode. Fig 3 illustrates the relationship between time and visibility set obtained by the 
circular array system at each angle. The horizontal axis stands for time and the vertical axis stands for rotating 
angle of the circular array system relative to its initial state. Each dot represents the visibility set, namely 
instantaneous samples obtained by the circular array system at a particular angle of the corresponding time. The 
blue dots indicate sample sets obtained by rotating circular array system. The red dots mean sample sets acquired 
by applying conjugate symmetry to samples represented by blue dots. As stated above, the circular array system 
needs to rotate half a circle to complete collecting all the visibilities needed for one image retrieval. A time 


sequence of visibilities for each sample set, or for each sample point, can be obtained after several imaging 
period, just as the horizontal solid line intersects 6 times with red and blue dots in fig 3 during the 6 imaging 
period, which means a time sequence of 6 values for each sample point in frequency domain. Then the 
interpolation method (nearest, linear, spline) can be used to get the interpolated value at a particular time from 
the time sequence of a sample point. Once the interpolation method is applied to each time sequence of all the 
sample points in frequency domain for a particular time t, then an image corresponding to time t can be retrieved 
from the interpolated samples at time t shown as vertical black dotted line in fig 3. 
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Fig 3. Relationship between time and visibility set obtained by the circular array system at each angle 


SIMULATION RESULT 


The interpolation method mentioned above is verified by GIMS simulation system, which can simulate 
time-sharing sampling process of GIMS, manipulate samples with the above interpolation method and retrieve 
the brightness temperature map with pseudo-polar IFFT retrieval algorithm [6]. The input to GIMS simulation 
system is a series of brightness temperature maps with 10s time interval generated by FNL/WRF/RTTOV 
method, which models brightness temperature variations of the full earth disk during the super typhoon phase of 
tropical cyclone (TC) Usagi. Fig 4(a) demonstrates the 50.3GHz input brightness temperature map at one time. 
The black ellipse illustrates the dynamic target—TC Usagi, which is zoomed in in fig 4(b). 


Given the 10 seconds spaced input brightness temperature map sequence and configured with imaging 
period of 5 minutes, the simulation system generates the blurred retrieved image corresponding to 5min imaging 
period and the de-blurred retrieved image by applying the interpolation method. Since Quality of retrieved image 
can be quantified using RMS error between retrieved image and input image, fig 5 shows RMSE value for 
tropical cyclone area of blurred image, de-blurred image and the snap shot image which amounts to retrieved 
image for static observational target and can be used as a contrast to illustrate impact of observing dynamic 
target. 


It can be discovered from fig 5 that RMSE of blurred image is higher than that of snap shot image for each 
of the 7 frequency channels, which means observing dynamic target indeed degrades the retrieved image. 
Comparing the three kind of interpolation method (nearest, spline, linear), nearest interpolation method makes 
no difference to the retrieved image. This might due to the fact that time sequence of samples for each sample 
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point is too sparse, making interpolated samples of nearest interpolation the same as the unprocessed samples of 
time-sharing sampling. However, spline and linear interpolation method can reduce the imaging error, especially 
the spline method, which has a relatively steady de-blurring effect for all the frequency channels as well as 
reducing image error to the level of the referential snap-shot image. 


(a) full earth disk {b) tropical cyclone Usagi area 
: 0.07 
TC Usagi: 

0.065 
g G 

2 2 0.06 
FA FA 
8 8 
5 5 

3 3 0.055 
A= E 
wv a an] 

0.05 

0.045 

-0.1 -0.05 0 0.05 0.1 0.15 0.02 0.025 0.03 0.035 0.04 0.045 
direction cosine (£) direction cosine (6 


100 120 140 160 180 200 220 240 260 280 


Fig 4. Input brightness temperature map at one time (50.3GHz, (a) full earth disk, (b) tropical cyclone Usagi 
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Fig 5. RMSE of snap-shot image, blurred image and de-blurred image for 7 frequency channels ((a) 50.3GHz, 
(b) 51.76GHz, (c) 52.8GHz, (d) 53.596GHz, (e) 54.4GHz, (f) 54.94GHz, (g) 55.5GHz) 


To illustrate effect of spline interpolation de-blurring method and exclude inherent imaging error in retrieval 
algorithm, fig 6 shows the difference between absolute error map of blurred retrieved image and that of 
de-blurred image after using spline interpolation method at 50.3GHz. It can be found spline interpolation reduce 
imaging error for tropical cyclone and the reduced error mainly distributes on the left part of tropical cyclone 
area, which coincides exactly with spiral rain band depicted in fig 4(b). This is easy to understand since the 
spiral rain band has more intense variation relative to the background when TC spins around its center, and thus 
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contributes most to the blur in the retrieved image. Therefore the de-blurring method will mainly have effect on 


this part. 
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Fig 6. Difference between absolute error map of blurred retrieved image and that of de-blurred image of 
spline interpolation method at 50.3GHz 


CONCLUSION 


The proposed interpolation method can reduce imaging error caused by time-sharing sampling of GIMS when 


observing dynamic target. Specifically, the spline interpolation method can de-blur the retrieved image to the 


level of referential snap-shot image, which corresponds to retrieved image for static observational target. 
Moreover, since the spline interpolation de-blurring method can be used to obtain restored image at any 
particular time, the method can also be employed to retrieve multiple images per imaging period, i.e. increase 


imaging temporal resolution. 
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